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Abstract 

In this paper, we introduce the Birnbaum-Saunders power series class of distributions which is ob- 
tained by compounding Birnbaum-Saunders and power series distributions. The new class of distributions 
has as a particular case the two-parameter Birnbaum-Saunders distribution. The hazard rate function 
of the proposed class can be increasing and upside-down bathtub shaped. We provide important math- 
ematical properties such as moments, order statistics, estimation of the parameters and inference for 
large sample. Three special cases of the new class are investigated with some details. We illustrate the 
usefulness of the new distributions by means of two applications to real data sets. 

Keywords: Birnbaum-Saunders distribution; Maximum likelihood estimation; Order statistic; Power 
series distribution. 



1. Introduction 



Birnbaum and Saunders (1969) introduced a two-parameter distribution in order to model the failure 
time distribution for fatigue failure caused under cyclic loading, called the Birnbaum-Saunders distribu- 
tion (BS). Since then, this distribution has been extensively used for modelling failure times of fatiguing 
materials in several fields such as environmental and medical sciences, engineering, demography, biolo- 
gical studies, actuarial, economics, finance and insurance. A random variable T following the BS(a,f3) 
distribution with shape parameter a > and scale parameter (3 > is defined by 
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where Z is a standard normal random variable. The cumulative distribution function (cdf) of T <~ 
BS(a,/3) has the form 

F BS (t;a,P) = $(V), t>0, (1) 

where v = a^ 1 p(t//3), p(z) — z 1 / 2 — z -1 / 2 and $(•) is the standard normal cumulative function. The 
corresponding probability density function (pdf ) is 



f BS {t; a, 13) = K (a, /?) t~ 3/2 (t + /?) exp 



r(t/r3) 
2a 2 



t > 0, 



(2) 



where K(a,j3) — cxp(a~ 2 )/ (2a \/2ir/3) and t(z) = z + z^ 1 . The parameter j3 is the median of the 
distribution, i.e. Fbs(P) — $(0) = 1/2- For any k > 0, kT <~ BS(a,kj3). Additionally, the reciprocal 
property holds, i.e., T~ x ~ BS(a, see Saunders (1974). The expected value and variance are 

E(T) = /3(1 + a 2 /2) and Var(T) = (a/3) 2 (l + 5a 2 /4), respectively. 

The BS distribution has a close relation to the normal distribution and has been widely studied in 
the statistical literature. For example, Ng et al. (2006) proposed a simple bias-reduction method to 
reduce the biases of the maximum likelihood estimators (MLEs) for the BS model. They also discussed 
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a Monte Carlo EM-algorithm for computing theses estimators. Kundu et al. (2008) investigated the 
shape of the BS hazard function. 

In the last few years, several classes of distributions were proposed by compounding some useful 
lifetime and power series distributions. Chahkandi and Ganjali (2009) studied the exponential power 
series (EVS) family of distributions, which contains as special cases the exponential geometric {EG), 
exponential Poisson {£ V) and exponential logarithmic (EC) distributions, introduced by Adamidis and 
Loukas (1998), Kus (2007) and Tahmasbi and Rezaei (2008), respectively. Morais and Barreto-Souza 
(2011) defined the Weibull power series (WPS) class of distributions which includes as sub-models the 
EVS distributions. The WVS distributions can have an increasing, decreasing and upside down bathtub 
failure rate function. The generalized exponential power series (QEVS) distributions were proposed by 
Mahmoudi and Jafari (2012) in a similar way to that found in Morais and Barreto-Souza (2011). In a 
very recent paper, Silva emphet al. (2012) introduced the extended Weibull power series (EWVS) class 
of distributions which includes as special models the EVS and WPS distributions. 

In a similar manner, we introduce the Birnbaum- Saunders power series (BSVS) class of distributions, 
which is defined by compounding the Birnbaum-Saunders and power series distributions. The main 
aim of this paper is to introduce a new class of distributions, which extends the Birnbaum-Saunders 
distribution, with the hope that the new distribution may have a 'better fit' in certain practical situations. 
Additionally, we will provide a comprehensive account of the mathematical properties of the proposed 
new class of distributions. The hazard rate function of the proposed class can be increasing and upside- 
down bathtub shaped. We are motivated to introduce the BSVS distributions because of the wide usage 
of ([I]) and the fact that the current generalization provides means of its continuous extension to still 
more complex situations. 

This paper is unfolds as follows. In Section 2, a new class of distributions is obtained by mixing 
the BS and zero truncated power series distributions. The mixing procedure was previously carried out 
by Morais and Barreto-Souza (2011) and Silva et al. (2012). In Section 3, some properties of the new 
class are discussed. In Section 4, the estimation of the model parameters is performed by the method 
of maximum likelihood. In Section 5, we introduce and study three special cases of the BSVS class 
of distributions. In Section 6, two illustrative applications based on real data are provided. Finally, 
concluding remarks are presented in Section 7. 



2. The new class 

Let Ti , . . . , T/v be independent and identically distributed (iid) BS random variables with cdf ([!]) 
andpdf {2]). We assume that the random variable N has a zero truncated power series distribution with 
probability mass function 

p n = P(N = n) = -^,n = l,2,..., (3) 

where a n depends only on n, C(9) = J2n=i a n and > is such that C(9) is finite. Table [lj reported 
in Morais and Barreto-Souza (2011), summarizes some power series distributions (truncated at zero) 
defined according to (J3| such as the Poisson, logarithmic, geometric and binomial distributions. 

By assuming that the random variables N and T are independent, we define X = min(Ti, . . . ,Tjv). 
Then, 

n f) n 

P(X<x,N = n) = -£r § j[l-(l-$(v)r], x>0, n>l. 
The BSVS class of distributions is defined by the marginal cdf of X: 

F BS vs{x; 9,a,p) = l- ° „~ ^ v))] , x > 0. (4) 

Hereafter, the random variable X following Q with parameters 9, a and f3 is denoted by X ~ 
BSVS{9,a,fi). The pdf to @ is 

f B svs(x; 0, a, ft = 9f BS (x) C , *>0, (5) 

where /bs(-) is given in 
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Distribution a n C(9) C'{6) C"{9) 0(e)- 1 6 

Poisson n!" 1 e e - 1 e e e° \og(9 + 1) £ (0, oo) 

Logarithmic n' 1 - log(l - (9) (1 - tf) -1 (1 - ff)- 2 1 - e" e (9 e (0, 1) 

Geometric 1 0(1 - f)" 1 (1 - £>)~ 2 2(1 - 9)~ 3 9(9 + l)" 1 6» 6 (0, 1) 

Binomial (™) (0 + l) m - 1 m(0 + l)™" 1 (7+l)2-i ^ " " 1 £ ^ 

Table 1: Useful quantities for some power series distributions. 

Proposition 1. The classical BS distribution with parameters a and f3 is a limiting special case of the 
proposed model when 8 — > 0. 

Proof. For x > 0, we have 

oo 

x>„ [0(1 -$(«))]" 

lim Fnsvs{%) = 1 — lim ^— — 

0^0+ 0^0+ 



n=l 

oo 

n=2 

oo 



= 1 - lim — = $(u). 

a 1 



n=2 



□ 

We now derive two useful expansions for the density function We have C'(9) — Y^n=i na n l 
and then 

f B svs(t;0,a,p) = ^(tja^J^n-^-Il-^u)]"- 1 . (6) 

n— 1 ^ ' 

By using the binomial theorem in equation ([6]), we can write 

oo n— 1 

fBSVs(x;0,a,f3) = f BS (x;a,f3)Y J ^2^n,k<l>(v) k , (7) 



71 = 1 fc = Q 



where u n ,k = {-!) np n 



n- 1 
k 



Proposition 2. The density function of X can be expressed as an infinite linear combination of densities 
of minimum order statistics of T . 

Proof. We know that C'(6) = J2n=i na n • Therefore, 

oo 

fBSVs(x;0,a,f3) = £p„ f T(1) (x;a, f3,n), x > 0, 

n=l 

where /rm (^j a j /3j ?i) is the density function of T^ = min(T 1 , . . . , T n ), for fixed n, given by 

/ T(1) (t; a, p,n) = n f BS (t)[l - $(«)]»-\ t > 0. 

□ 
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Remark 1. Let Y — max(Ti, . . . , TV), then the cdf and pdf of Y are 

F Y (y;6,a,P) = ^^- , y>0 (8) 

and 

f Y {y, 6, a, p) = 6f BS (y) c ' { ^ )} , y>0. 

Note that 

OO 

f Y (y,0,a,f3) = ^2p n f T(n) {y;a,l3,n), y > 0, 

n=l 

where /t (?i) (y; c 5 Pi n) is the density function of T( n ) = max(T 1; . . . , T n ), for fixed n, given by 

f T{n) (t; a, /?, n) = n / 85 (t) S^)"" 1 , i > 0. 

The distribution with cdf ([8]) is called the complementary Birnbaum- Saunders power series (CBSVS) 
distribution. This class of distributions is a suitable model in a complementary risk problem based in the 
presence of latent risks which arise in several areas such as public health, actuarial science, biomedical 
studies, demography and industrial reliability (Basu and Klein, 1982). However, in this work, we do not 
focus on this alternative class of distributions. 

The BSVS survival function becomes 

W*;W) = m^M , x>0 

and the corresponding hazard rate function reduces to 

h BSVS (x; 6, a, p) = 9f BS (x) C ffi 1 ~ f ( " ))] , x > 0. 

C[0{1 - 9{v))\ 

3. Quantiles, order statistics and moments 

The quantile function of X, say Q(u), is given by 



! = (){„)= + 
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where 5 = — $ 1 {C 1 [(1 — u)C(9)} /#}, 4> : (-) is the inverse cumulative function of the standard normal 
distribution, C _1 (-) is the inverse function of C(-) and u is a uniform random number on the unit interval 
(0,1). 

Let Xx , . . • , X m be a random sample with density function ([5| and define the ith order statistic by 
T im . The density function of T i:m , say f i:m (x;6,a, p), is given by 

f i:m (x;8,a,P) = — ] ——zfBSVs(x;0,a,P)F BSVS (x;9,a,py^ 1 [1 - F BSVS (x; 9, a, P)] m ~ l 

B(i, m — i + 1) 



f BS {x',a,0) sr^f .aA-A, , a 



-B(z, m — i + 1) 



fc=0 



where Sgsvsi') is the BSVS survival function. 

Proposition 3. The BSVS order statistics can be expressed as a linear combination of exponentiated 
Birnbaum- Saunders (EBS ) distributions with parameters a, P and s + 1. 

In fact, 

i— 1 oo m+k—i+j 

fi- m {x]9,a,P) =X]X] X! w fc,J>' l s+l( a; > Q! >/ ? )) 
fc=0 3=0 s=0 



where 



(s + l)B(i, m-i+ l)C*(6»)™+ fe - i ' 
1 r 

do = a" l+k ~\ d r = — 1~ E [ J '( m + k-i)-r + j] a j+ id r -j 



rd . 
3=1 



and h s +i(t, a, P) denotes the EBS density function with parameters a, ft and s + 1. 

Many of the important characteristics and features of a distribution are obtained through the mo- 
ments. The ordinary moments of X can be derived from the probability weighted moments (PWMs) 
(Greenwood et al, 1979) of the BS distribution formally defined for p and r non- negative integers by 



K(a,P) / i p ~ 3/2 (i + ^)exp 



r(t/P) 
2a 2 



<$> r (v)dt. (9) 



The integral M can be easily computed numerically in software such as MAPLE, MATHEMATICA, 
Ox and R, see Cordeiro and Lemonte (2011). An alternative expression for r Pjr . can be derived from 
Cordeiro and Lemonte (2011) as 

w = ?E(;) E A (fcl ,...,^) 2 E(-ir(\ +j 

3=0 V/ ki,...,kj=0 m=0 V 

x pr ( 2 «i+3-2m)/2 j( p + ( 2g , +j _ 2m )/2, a) . 

Here, Sj = fc x + • • • + A(fc X) . . . , k/) = a~ 2s ^a kl . . . a kj , a fc = {-l) k 2 ( - 1 - 2k: ^ 2 {^{2k + l)^.}- 1 
and 

w . _ K p+1/2 (g- 2 ) + if p -i/ 2 (a:~ 2 ) 
1 (p, a) — , 

where the function K v (z) denotes the modified Bessel function of the third kind with v representing its 
order and z the argument. Its integral representation is K v (z) — 0.5 J_ exp{— z cosh(t) — i>t}dt. 
Thus, the sth moment of X can be obtained from equations Q and ^ as 

oo n— 1 

E{X S ) = E£>„, fe T s>fe . 

n=l fc=o 

Expressions for the sth moments of the order statistics Xi :m , . . . , X m>m with cumulative function Q 
can be obtained using a result due to Barakat and Abdelkader (2004). We have 

E(X? m ) = s jr (_l)i-m+i-l ^ ^ ^ ^ ^ 



DC- 



j=m— «+l 

for i = 1, . . . , to. 

4. Maximum likelihood estimation 

Here, we discuss maximum likelihood estimation and inference for the BSVS distribution. Let 
Xi, . . . ,x n be a random sample from A and let = (9, a, (3) T be the vector of the model parameters. 
The log-likelihood function for reduces to 

o " n 

£(9, a, p) = n{\og{6) - log[C(0)] + log[ K (a, P)}} - - £ logfc) + £ \og( Xl + p) 

i=l i=l 



1 n / \ " 

3^ 7 [*(!-*(«,))]}. 



2a 2 ^ \P 

i=i v K 
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The components of the score vector U (0) are given by 



n C'{6) , ^ l-<f>(vi) dC [0(1 - $( Vi ))] 



n _ C"(9) 1- 

e{ ~'~ 6 n C(9) + f-^ & [9(1 - 90 

Me) - = f i + A) + A ±r ( 5) + 4 ± 3C ' - •"»»' 

v ; a\ a 2 / a 3 f-' V /3 / a 3 ^ C" 6>(1 



x /3 7 a 3 ^ C [0(1 - $( Wi ))] 9a 



„ , , x /3\ * <KVi)T(yfc/p) dC'[9(l 



and 



2/3 2a 2 /3^'V/ 9 2a/3 -f^ C" [0(1 - $(i>i))] 9/3 

where </>(•) is the standard normal density and u< = a~ l p(xi/9). 

Setting these equations to zero, U(&) = 0, and solving them simultaneously yields the MLE of 
0. These equations cannot be solved analytically and statistical software can be used to compute them 
numerically using iterative techniques such as the Newton-Raphson algorithm. 

Often with lifetime data and reliability studies, one encounters censoring. A very simple random 
censoring mechanism that is often realistic is one in which each individual i is assumed to have a lifetime 
Xi and a censoring time Ci, where Xj and Ci are independent random variables. Suppose that the data 
consist of n independent observations Xi — min(X;, Ci) and 5i = I(Xi < Ci) is such that 8i = 1 if Xi is 
a time to event and 5, = if it is right censored for i = 1, . . . , n. The censored likelihood L(9, a, /3) for 
the model parameters is 



L(9, a, /3) oc Y[ [fBSVs(xt;9, a, f3)] Si [S BSVS (xi; 9, a, /?)] 



l-Si 



In order to perform interval estimation and hypothesis tests on the model parameters a, /3 and 9, 
the normal approximation for the MLE can be applied. Note that under certain conditions for the 
parameters, the asymptotic distribution of y/n(& — 0) is multivariate normal A^(0, 1~ (©)), where 
J n (0) is the observed information matrix and I(&) = linin^oo J„(0). Note also that an 100(1 — 7)% 
(0 < 7 < 1/2) asymptotic confidence interval for the ith parameter Bi in is specified by 



ACh = (Bi - Zl _ 7/2 Vj^, Bi + Zl _ 7/2 y/fc 



where J 9i > Bi stands for the ith diagonal element of the inverse of the observed information matrix esti- 
mated at 0, i.e., J„(0) _1 , for i = 1,2,3, and z 1 f / 2 is t ne 1 — 7/2 standard normal quantile. 

For testing the goodness-of-fit of the BSVS distribution and for comparing it with some sub-models, 
we can use the likelihood ratio (LR) statistics. The LR statistic for testing the null hypothesis T-Lq: 
©i = 0^ versus the alternative hypothesis Hi: ©i ^ ©^ is given by w — 2{£(&) — £(©)}, where 
and are the MLEs under the alternative and null hypotheses, respectively. Under the null hypothesis, 
Ho, ro-4x^, where k is the dimension of the subset ©i of interest. 



5. Special cases 

In this section, we investigate some special models of the BSVS class of distributions. We offer some 
expressions for the moments and moments of the order statistics. We illustrate the flexibility of these 
distributions and provide plots of the density and hazard rate functions for selected parameter values. 

5.1. Birnbaum-Saunders geometric distribution 

The Birnbaum-Saunders geometric (BSQ) distribution is defined by the cdf Q with C(9) — 9(1 — 9)~ 1 
leading to 

F BS g(x; 6, a, p) = 1 - ( * ~ gjj^S > * > 0, (10) 

where 8 G (0,1). 



Proposition 4. The distribution of the form (10) is geometric minimum stable. 

Proof. The proof follows easily of the arguments given by Marshall and Olkin (1997, p. 647). We omit 
the details. □ 

The associated density and hazard rate functions are 

f BSG (x; 0, a, 13) = (1 - 6)f BS {x) [1 - 9(1 - ^(v))}' 2 

and 

h ( a R \ Ibs{x) h BS (x) 

kB8 O {xA«,0) = [!_*(„)]{!_,,[! = _*(„)]' ^ 

for x > 0, respectively. From equation (111, we note that h B sg(x;9,a, (3) /h B s(x;a, j3) is decreasing in 
x. It can also be shown that 

1 



lim hBsg(x\6, a,/3) = and lim h B sg(x; 9, a, /?) . . 

a;— >0 a— >oo let p 

Corollary 1. Lei X - BSQ(9,a,f3) and T ~ BS{a,f3). Then h Bsg (x; 9, a, (3) > h BS (t;a,(3) for all 

0e(O,l). 

Proof. It is straightforward. □ 

Cancho et al. (2012) proposed and studied the geometric Birnbaum-Saunders regression with cure 
rate. However, they considered C(9) = (1 — 9)~ 1 instead of C{9) = 9(1 — 9)~ x . In Figure[lJ wc plot the 
density and hazard rate functions of the BSQ distribution for selected parameter values. We can verify 
that this distribution has an upside-down bathtub or an increasing failure rate function depending on 
the values of its parameters. 

Expressions for the density function and the moments of the order statistics Xi :m , . . . , X m:m from a 
random sample of the BSQ distribution are given by 



i—X /. -, \ r /-i n\ ^ / \-im-\-k—i 

(1 - 9)&(-v) 1 



1 - 6$(-v) 



±s(i,m — i + l ^-^ V k 
v ' ; fe=0 v 

and 

which are easily determined numerically. 

Let x\, • • • , x n be a random sample of size n from X ~ BSQ (8, a, /?). The log- likelihood function for 
the vector of parameters = (9, a, f3) T can be expressed as 

£(&) = n{log(l-0)+logK«,/3)]}--^log(x l )+^log(x i +/3) 

i=i i=i 

^EMl )-2£]og[l-*(l -*(«,))]. 
i=l v ^ / i=1 

5.2. Birnbaum-Saunders Poisson distribution 

The Birnbaum-Saunders Poisson (BSV) distribution follows by taking C(9) = e e — 1 in Q, which 
yields 

jr a a mi exp{fl(i - 3>(u))} - l 

FBSv{t;9,a,f3) = 1 — . 

exp(0J — 1 

The density and hazard rate functions of the BSV distribution are 

hsv(t,9,a,[3) = eM9) _ 1 

and 

, nn^ 9f BS (t)eM0(l-Hv))} 
h BSV (t,9,a,P) - exp{0(1 _ $(u))} _ 1 ■ 

7 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 
f 



(a) a = 0.2 and /3 = 5.0 



(b) a = 2.0 and /3 = 5.0 
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Remark 2. The limit of the BSV hazard rate junction as t — > is 0. 



In Figure [2] we plot the density and hazard rate functions of the BSV distribution for selected 
parameter values. This distribution can have an upside-down bathtub or an increasing failure rate 
function depending on the values of its parameters. 




Figure 2: Plots of the BSV density and hazard rate functions for some parameter values. 



The expressions for the density function and the moments of the order statistics Xi :m , . . . ,X m:r 
from a random sample of the BSV distribution are 

fBs(x;a,l3) ^ k fi - l\ J expjg^)^-! 1 m+fc -' 
/,:ml J "B(i,m-i + l)^ ' \ k )\ e e -l j 

and 

™ C_i y'-m+i-l / ■ _ i \ / \ roo 
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Let xi t ■ ■ ■ , x n be a random sample of size n from X ~ BS'P(9, a, f3). The log-likelihood function for 
the vector of parameters — (6, a, f3) T can be expressed as 

q n n 

£(@) = n{\og(0) - log[exp(0) - 1] + log[«(a, (3)}} - - g logfo) + £ log(s, + /3) 



2 

t=l »=1 

i=l v y 7 i=l 

5.5. Bimbaum- Saunders logarithmic distribution 

The cdf of the Birnbaum-Saunders logarithmic (BSC) distribution is defined by Q with C(9) 
— log(l — 0), < < 1, which corresponds to the logarithmic distribution. We obtain 



The associated density and hazard rate functions are 



fssc(x;O,a,0) 



Iog(l-0){1-0[1 -#(«)]} 
and 

Remark 3. TTie Zzrmi o/ £/ie BSC hazard rate function as x — ¥ is 0. 

In Figure [3j we plot the density and hazard rate functions of the BSC distribution for selected 
parameter values. We can verify that this distribution can have an upside-down bathtub or an increasing 
failure rate function depending on the values of its parameters. 

The expressions for the density function and the moments of the order statistics Xi :m , . . . ,X m:m 
from a random sample of the BSC distribution are 

h s (x;a,/3) ^ fe (i l\ f log [1 - 8<£(-v)} \ m+fc -' 
U -- m[X) - B (i,m-i + l)^ [ L) \ k J\ log(l-0) J 

and 



j—m—i 



+1 [iogC 1 - W \ m -y \.i J Jo 



1^3 



1 - 0$(-t>) 



dx. 



Let cci, • • • , i„ be a random sample of size n from X ~ BSC(9, a. f3). The log-likelihood function for 
the vector of parameters = (8, a, /3) T can be expressed as 

£(&) = n {log(0) - log [- log(l - 0)] + log[«(a, /?)]} - - £ M^) + E lo S^ + 0) 



2 

i=l i=l 



1 " / \ - 



2a 2 ^ V 



6. Applications 

In this section, we compare the fits of some special models of the BSVS class by means of two real 
data sets to show the potentiality of the new class. In order to estimate the parameters of these special 
models, we adopt the maximum likelihood method (as discussed in Section 4) and all the computations 
were done using the subroutine NLMixed of the SAS software. Obviously, due to the genesis of the BS 
distribution, the fatigue processes are by excellence ideally modeled by this distribution. Thus, the use 
of the BS and BSVS distributions for fitting these two data sets is well justified. 
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Figure 3: Plots of the BSC density and hazard rate functions for some parameter values. 
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First, we consider a data set from Murthy et al. (2004) consisting of the failure times of 20 me- 
chanical components. The data are: 0.067, 0.068, 0.076, 0.081, 0.084, 0.085, 0.085, 0.086, 0.089, 0.098, 
0.098, 0.114, 0.114, 0.115, 0.121, 0.125, 0.131, 0.149, 0.160, 0.485. The MLEs of the parameters (with 
corresponding standard errors in parentheses), the value of —2^(0), the Kolmogorov-Smirnov statistic, 
Akaike information criterion (AIC) and Bayesian information criterion (BIC) for the BSG, BSV, BSC 
and BS models are listed in Table [2) Since the values of the AIC and BIC are smaller for the BSQ 
distribution compared with those values of the other models, the new distribution seems to be a very 
competitive model for these data. 

Table 2: Parameter estimates, K-S statistics, AIC and BIC for failure times data 



Distribution 


a 





6 


K-S 


-2£(@) 


AIC 


BIC 


BSQ 


0.6461 


0.4521 


0.9950 


0.1314 


-77.6 


-71.6 


-68.6 




(0.6194) a 


(0.8247) 


(0.0184) 










BSV 


0.4774 


0.1735 


5.1057 


0.1224 


-73.2 


-67.2 


-64.3 




(0.0877) a 


(0.0304) 


(2.0932) 










BSC 


0.3549 


0.2437 


0.9999 


0.2055 


-74.0 


-68.0 


-65.0 




(0.0570) a 


(0.0419) 


(0.0001) 










BS 


0.4466 


0.1107 




0.2029 


-65.5 


-61.5 


-59.5 




(0.0706) a 


(0.0108) 













a Denotes the standard deviation of the MLE's of a,/3 and 9. 

Plots of the pdf and cdf of the BSQ, BSV, BSC and BS fitted models to these data are displayed 
in Figure |4j They indicate that the BSG distribution is superior to the other distributions in terms of 
model fitting. Additionally, it is evident that the BS distribution presents the worst fit to the current 
data and then the proposed models outperform this distribution. 





(a) 



(b) 



Figure 4: Estimated (a) pdf and (b) cdf for the BSQ, BSV, BSC and BS models for failure times data. 



We can also perform formal goodness-of-fit tests in order to verify which distribution fits better to 
these data. Chen and Balakrishnan (1995) proposed a general approximate goodness-of-fit test for the 
hypothesis Hq : X\, . . . ,X n with Xi following F(x;9), i.e., under Hq, Xi, . . . , X n is a random sample 
from a continuous distribution with cumulative distribution F(x;6), where the form of F is known but 
the p-vector 6 is unknown. The method is based on the Cramer- von Mises (C-M) and Anderson-Darling 
(A-D) statistics and, in general, the smaller the values of these statistics, the better the fit. 
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Table[3] gives the values of the C-M and A-D statistics (p- values between parentheses) for the failure 
times data. Thus, according to these formal tests, the BSQ model fits the current data better than the 
other models. These results illustrate the potentiality of the BSQ distribution and the importance of the 
additional parameter. 



Tabic 3: Goodncss-of-fit tests 



C-M A-D 



Bsg 


0.0469 (0.5563) a 


0.3116 (0.5517) a 


BSV 


0.0877 (0.1650) 


0.6629 (0.0835) 


BSC 


0.0784 (0.2173) 


0.6064 (0.1151) 


BS 


0.1967 (0.0059) 


1.3748 (0.0015) 


a Denotes the p- value of the test. 



As a second application, we shall analyze a data set from McCool (1974) on the fatigue lives (in 
hours) of 10 bearings of a certain kind. The data are: 152.7, 172.0, 172.5, 173.3, 193.0, 204.7, 216.5, 
234.9, 262.6, 422.6. Table [4] gives the MLEs (standard errors between parentheses) of the parameters, 
the values of —2^(0) and of the Kolmogorov-Smirnov, AIC and BIC statistics for the fitted BSQ, BSV 
and BS models to the second data set. Plots of the fitted BSQ, BSV and BS densities are given in 
Figure [5] 



Table 4: Parameter estimates, K-S statistics, AIC and BIC for fatigue life data 



Distribution 


a 


d 


e 


K-S 


-2£(®) 


AIC 


BIC 


BSG 


0.3087 


350.98 


0.9672 


0.1681 


106.9 


112.9 


113.8 




(0.1285) a 


(182.50) 


(0.0861) 










BSV 


0.2917 


259.20 


3.1140 


0.1633 


108.3 


114.3 


115.2 




(0.0772) a 


(44.4148) 


(2.4589) 










BS 


0.2825 


212.05 




0.1707 


109.9 


113.9 


114.5 




(0.0632) a 


(18.7530) 













a Denotes the standard deviation of the MLE's of a, ft and 6. 
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The figures in Table [4] indicate that the BSQ model yields a better fit to these data than the other 
models. Table [5] gives the values of the C-M and A-D statistics (p- values between parentheses) for the 
fatigue life data. So, these values indicate that the null hypothesis is strongly not rejected for the BSQ 
distribution, whereas the null hypothesis is weakly not rejected for the other models. Thus, according 
to these goodness-of-fit tests, the BSQ model fits the current data better than the other models. 



Tabic 5: Goodncss-of-fit tests 



Model 



Statistics 



C-M A-D 



BSQ 0.0370 (0.7373) a 0.2761 (0.6575) 3 

BST 0.0573 (0.4108) 0.4243 (0.3178) 

BS 0.0862 (0.1725) 0.6148 (0.1098) 



Denotes the p- value of the test. 



7. Concluding remarks 

The two-parameter Birnbaum-Saunders distribution is extended to a new class of fatigue life dis- 
tributions by compounding the Birnbaum-Saunders and power series distributions, thus creating the 
Birnbaum-Saunders power series distribution. As expected, the new class includes as special model the 
BS distribution. The BSVS failure rate function can be increasing and upside-down bathtub shaped. 
Some mathematical properties of the new class are provided such as the ordinary moments, quantile 
function, order statistics and their moments. The estimation of the model parameters is approached by 
the method of maximum likelihood and the observed information matrix is derived. We fit some BSVS 
distributions to two real data sets to show the potentiality of the proposed class. 

Future research could be adressed to study the complementary Birnbaum-Saunders power series 
distribution. This class of distributions is a suitable model in a complementary risk problem based in 
the presence of latent risks which arise in several areas such as public health, actuarial science, biomedical 
studies, demography and industrial reliability. 
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Appendix A. Elements of the observed information matrix 

The observed information matrix </(©) for the parameters 9, a and /3 is given by 

F) 2 P(Pt\ I U 09 ^ 9a ^8$ \ 
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whose elements are 
Uee 



U 6c 



UryOL 



n 



n 



C"{9) _ ( C'(0) Y 

c(e) {c(e)J 



a 



C' a {-0*{-v)) 
C>(-6$(-v)) ' a 2 ^ 



c"(-o$(-vi)y 

Of/ 

C r (-r9$(-w)) 



1 



2a/3 C (-6<f>(-v)) 2a/3^ C" (-0$(-t>)) 2af3 j-^ 



C"(-0$(-v)) 

c(-e&(-v)) ' 



a* 



1 + 



2g " c; (-^(-u)) 

3 /7' ' 



Q: 



~ > 4>{vi)p - C' (-0$(-u)) — 

1 ™ 



rf C" ' a 2 ^ C'(-0#(-u)) 



2 ' 



6 



)) 



-jc Q (-(?$(-«)) C / M$( _ w)) 



+ 





a 2 ^ C" (-(9$(-v)) : 



1 



2/3 2 2a 2 £ 2 

g 2 
2a 2 /3 



n 



E 



=i 

2_ 



+ ^ 

*i 



2a 2 /3 ; 



E- 

»=i 



2^ ry 



i=i 



+ 



+ 



i=l 



fU\ ( U_ fl_ ^ FtA C fi (-m-v)) 
Pj\2pyti 2t t \jp) O (-6<S>(-v)f 

fU , [6\C' fi (-0*(-v)) 



Ul C (-0$(-v)) 4a/3 2 ^ 



E^^ 



2a/3 ^ 



E^(^)' 1 



i=i 



t"\ C' (-*$(-,;)) 



/3 / C"(-6»$(-w)) ' 



i"\ C' (-*<!>(-«)) 



/3 / (7 (-0$(-v)) 



where C^(-) is the first derivative of C'(-) with respect to fc and Cj fe (-) corresponds to second derivative 
with respect to j and fc. 
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